Manifestation of the Arnol'd Diffusion in Quantum Systems 
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We study an analog of the classical Arnol'd diffusion in a quantum system of two coupled non- 
linear oscillators one of which is governed by an external periodic force with two frequencies. In 
the classical model this very weak diffusion happens in a narrow stochastic layer along the coupling 
resonance, and leads to an increase of total energy of the system. We show that the quantum 
dynamics of wave packets mimics, up to some extent, global properties of the classical Arnol'd 
diffusion. This specific diffusion represents a new type of quantum dynamics, and may be observed, 
for example, in 2D semiconductor structures (quantum billiards) perturbed by time-periodic external 
fields. 
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As is known, the main mechanism for the onset of dy- 
namical chaos in classical Hamiltonian systems is the de- 
struction of separatricies of non-linear resonances due to 
perturbation terms. However, for a small enough per- 
turbation in conservative systems with two degrees of 
freedom (N = 2), as well as for one-dimensional time- 
dependent Hamiltonians, chaotic regions are bounded 
in the phase space by the Kolmogorov-ArnoPd-Moser 
(KAM) surfaces Q. In this case the chaos is non-global 
since chaotic trajectories are located in restricted regions 
of the phase space of a system. 

The situation changes drastically for many- 
dimensional systems (N > 2), where KAM surfaces 
no longer separate one stochastic region from another, 
and chaotic layers generically form a stochastic web [Q. 
Thus, if trajectory starts in the vicinity of a specific non- 
linear resonance, it can diffuse along stochastic layers 
of many resonances that cover the whole phase space. 
This universal global instability known as the Arnol'd 
diffusion ^ , for the first time was observed in numerical 
experiments J 4[, and then studied in detail in (see 
also review j7] and the book 0). 

It is very difficult to observe this specific diffusion since 
its rate is exponentially small, and it manifests itself only 
for initial conditions inside very narrow stochastic layers. 
However, it is believed that the Arnol'd diffusion may 
play an important role in different physical situations. 
For example, it is argued that this diffusion is crucial for 
the estimate of the stability of our Solar system, and is 
responsible for a loss of electrons in magnetic traps (see 
in [7]]). Also, the Arnol'd diffusion may have a strong in- 
fluence for the dynamics of protons in high energy stor- 
age rings. Recently, the possibility of observation of the 
Arnol'd web for a Hydrogen atom in crossed electric and 
magnetic fields has been discussed in Ref . || . 

So far, all the studies of the Arnol'd diffusion refer to 
classical models. On the other hand, it is important to 
understand what is the fingerprint, if any, of this dif- 
fusion in quantum systems. This question is far from 
trivial since quantum effects are known to typically sup- 
press classical effects of an exponentially weak diffusion 



H . In this Letter we perform a detail study of the quan- 
tum Arnol'd diffusion in a simple model, and show what 
are peculiarities of quantum dynamics that are due to 
this phenomenon. 

The system under consideration consists of two cou- 
pled quartic oscillators one of which is perturbed by a 
two-harmonic force. Following to Refs. || where the cor- 
responding classical model was studied both analytically 
and numerically, wc write the Hamiltonian, 



H = H° + H° - (ixy - f x{cos £2i* + cos Cl 2 t). (1) 
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and /j,, fo, fti, f^2 stand for the coupling constant, am- 
plitude and two frequencies of an external force, respec- 
tively. For momentum and coordinate operators the stan- 
dard commutation relations are assumed, [p x , x] = —ih~o, 
\p v , y] = —iho, with Hq as the dimensionless Plank con- 
stant. 

In the classical model the separatrix of the main cou- 
pling resonance oj® ~ u>y (determined in the absence of 
perturbation, fo = 0) is destroyed and creates a very 
narrow stochastic layer. The Arnol'd diffusion is caused 
by two driving terms with commensurate frequencies fii 
and £1%, that "force" a chaotic trajectory diffuse along 
the resonance layer . 

To consider the quantum model, first, one should find 
quantum eigenstates corresponding to this coupling res- 
onance, for fo = 0. Then, by switching on the perturba- 
tion, the analysis can be performed by using the Floquct 
formalism, since the perturbation is periodic in time. In 
this way, one can construct the evolution operator that 
allows one to study the dynamics of the model. 

In order to find stationary states corresponding to the 
coupling resonance, we write the wave function in the 
form, 



ip(x, y) 



(3) 
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Here ipn( x ) an d i^miv) are the eigenfunctions of H® and 
Hy (calculated numerically), and coefficients c niTn satisfy 
to the stationary Schrodinger equation, 

Ec n)7r i (E n -\- E ln )c nm fl ^ ^ %n,n f ym,m' Cn f ,m' ■> (4) 



where E n and E m are the eigenvalues of the Hamiltoni- 
ans H° and H°. 

y 

Our interest is in the dynamics of the model in the 
vicinity of the main coupling resonance determined by 
the condition uj no = uj„ 1q . Here h uj no = E' no , h Lo mo = 
E' mol and no = mo defines the resonance centre. In this 
region one can expand E n and E m in the Tailor series, 
keeping second order terms. It is convenient to intro- 
duce the indexes p = k + I and k, where k = n — no and 
I — m — tuq. Then, one can write, 



Ec Kp = h up + E^ Q (k 2 -pk+£ 



<•(••• + J2k' x k,k>y P -k-l-k' Ck'-V 
+ J2k> Xk : k'Vp-k-k' Cfc',o + 
+ J2k' x k,k'Up-k,l-k' Ck',1 + • 



(5) 



with lo = LJ na . Since matrix elements 

Xm,n and ym,n 

of coordinates are equal to zero for transitions between 
states of the same parity, the solution of Eqs. (||) consists 
of two independent sets, for odd and even values of p. 
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FIG. 1. Energy spectrum of the system ^ in units of 
Tiqlu. Five groups with 121 states in each group are shown 
for n = 10" 4 , h « 1.77 ■ 10~ 5 , and n = 446. 



uj is the classical frequency of small phase oscillations at 
coupling resonance. Accumulation points correspond to 
classical separatricies, and all energy levels inside sepa- 
ratricies are non-degenerate. The states slightly above 
separatricies are quasi-degenerate due to the symmetry 
of a rotation in opposite directions. 

In accordance with the spectrum structure it is conve- 
nient to characterize the states at coupling resonance by 
two indexes, by the group number q and by the index s 
which stands for the levels inside groups. Then, the en- 
ergy in each group can be written as E g<s = Ti^Loq + E qg , 
where E¥ s is the Mathieu-like spectrum of one group. 

Now we analyze the dynamics of the model in the pres- 
ence of the external force. For commensurate frequencies 
Q\ and Q2, the perturbation is periodic with the period 
T = iTx = jT 2 , where 7\ = 2tt/Oi, T 2 = 2tt/J1 2 and 
i, j are integers. The initial conditions were taken for 
a system to be about half-way between the two driving 
resonances, u = (fli + SI2) /2. 

Since the Hamiltonian ([I]) is periodic in time, the so- 
lution of the Schrodinger equation can be written as 



~>P(%,y,t) = exp 



' h J 



UQ{x,y,t). 



(6) 



Here UQ(x,y,t) = UQ(x,y,t + T) and eq are the 
quasienergy (QE) functions and quasienergies, respec- 
tively. They are determined by the evolution operator 
U describing the dynamics of our system in one period 
of the external field, 



UuQ(x,y) = exp 



ie Q T 



UQ{x,y). 



(7) 



Here the argument t is omitted since we are interested 
in the wave function only for discrete times NT, with 
N = 1,2,.... 

It is now naturally to represent the QE functions in 
the form u Q (x,y) = Y^ q , s A q, s ^qA x ^ v)> where ip q , s (x,y) 
are eigenstates of the unperturbed Hamiltonian H = 
+ Hy — [ixy. The coefficients A® s are the eigen- 
vectors of the operator U in the representation of two 
coupled nonlinear oscillators, that can be found by diag- 
onalization of the corresponding matrix U q ^ q '. s i. This 
matrix can be numerically obtained in the following 
way. Let the evolution operator U act on the initial 

(90, s ) 



Let us now consider the case of a small nonlinearity, state C q J' (0) — 5 q ^ qg 5 s , Sg . Then the wave function 
when the condition h Q cup > E% [k 2 -kp+^j is satis- C q q °' So) (T) at time T forms the column of the evolution 

fled. The corresponding numerical data for a fragment of 
the energy spectrum is shown in Fig.[l]. One can see that 
the spectrum consists of series of energy levels, that are 
shifted one from another by the value h u). The structure 
of energy spectrum in each group is typical for a quan- 
tum nonlinear resonance Go]. Lowest levels are practi- 
cally equidistant with the spacing equal to Houj, where C, 



(90, so 
q,s 

operator matrix, 

U q , s , ql ,AT)C^f o) (0) = U 9 ,. m ,. (T) = C q T o) (T). 



(8) 



Repetition of this procedure for different initial states 



(<zV) 

q,s 



(0) 



5q,q'5 s 



determines the whole matrix 
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U qjS;ql ^(T). As a result, the wave function C$' so) (T) 
can be computed numerically by integration of the non- 
stationary Schrodinger equation, 



ih C q , s = (houq + E*{) C a . R - 



q>,s' •"?,»)?'■«' 



q,sj ^q-s 
(COS flit + COS Q,2t) C q i <s i 



(9) 



If we introduce slow amplitude b qjS (t) via the transfor- 
mation 



C q , s (t) = b q , s (t) exp [-t (quj + E^/ho) t] , 
then after some algebra one can obtain 
ift b g<s = -/ocos (^t) x 

~V \ U -i(E M ,, ,-E^ B )t/h 



(10) 



(11) 
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where <5f2 = fix — £1%. Using the resonance approxima- 
tion, we keep only the most important slowly oscillating 
terms with q' = q ± 1 . 

The matrix elements x q , s - q ±i iS > in ( |TT| ) define transition 
probabilities between the states s and s' from the neigh- 
bor groups. This matrix has a specific block structure 
shown in Fig.Ha. 
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FIG. 2. (a) Matrix elements |Xg, s; g±i ]S '| that define the 
transition probability alongthe coupling resonance, for the 
same parameters as in FigM (b) Detailed structure of the 
"cross" -like broadening. 

Black and white regions correspond to big and small 
matrix elements, respectively. One fragment of this 
structure is shown in Fig.|]b in more details. The matrix 
elements in the center of Fig.||b correspond to the tran- 
sitions between lowest states (at the center of coupling 
resonance) of the groups with q — and q = 1. Matrix 
elements at the corners of Fig.^Jb, where s, s' > 30 or 
s, s' < —30, define transitions between the states above 
separatrix. All these elements quickly decrease with an 
increase of the difference \s — s'\. The "cross" -like broad- 
ening in Fig.^Jb, where matrix elements are large, corre- 
sponds to the transitions between separatrix states. As 
a result, the transition probability between separatrix 



states of neighbor groups (along the coupling resonance) 
is much larger than the transition probability between 
other states. This phenomenon is analogous to the quan- 
tum diffusion inside a separatrix, which was investigated 
in a degenerate Hamiltonian system pT| . 

Solving numerically Eqs.(p"l|), we obtain the matrix 
Uq, S ;q' ,s'{T) that determines the eigenvalues eq and 
QE functions A® s . Global properties of QE func- 
tions can be understood in terms of their "centers" 
E 9E S \ A q,s\ 2 , and dispersion a q = £„(<?- 



|^4? s | 2 in the unperturbed (fo = 0) basis, see 



Q = 

Fig.y. One can see that for a relatively small fj, and fo, the 
QE functions do not couple unperturbed states with dif- 
ferent values of q. On the other hand, a stronger coupling- 
results in a kind of "derealization" that is characterized 
by the spread of QE functions over many groups of states 
with different q. This fact manifests a large probability 
for the transition between those states that are involved 
in the diffusion along the coupling resonance. 
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FIG. 3. Plot of q versus y[Sq~ for QE functions A® s in 
the regions q — 0,±1 for different coupling constants: a) 
H = 3 • 10~ 5 , b) fx = 10~ 4 . Points correspond to different 
QE functions. 



Direct computation of the evolution of wave packets 
has confirmed the occurrence of a weak diffusion in en- 
ergy space, that is analogous to the classical Arnol'd dif- 
fusion. Numerical data have been obtained by computing 
the evolution matrix according to the relations, 



U a 



,(NT)=J2 A Z A q>'s>™P 



NT 



h 



CqANT) 



2_j Uq,s;q',s' 
q',s> 



(JVT)CV, fl ,(0). 



(12) 



(13) 



We have studied quantum dynamics and calculated dif- 
fusion coefficient for different initial states in the regime 
when the values of fi and fo are small enough and 
/o//i = 0.01, so that the coupling and two driving reso- 
nances do not overlap. 

Quantum dynamics for different initial conditions is 
illustrated in Fig.^ where the energy dispersion A q = 
(Al?)7^ 2 with A q = V (g - qfZs 



\C a 



and 



3 



q = J2q1^2s |Cg,s | 2 is plotted versus the rescaled time 
N = t/T. The figure clearly shows different character 
of the evolution for three initial states. For the states 
taken from below and above the separatrix, the energy 
width of packets oscillates, in contrast with a diffusion- 
like time-dependence for the case when the initial state 
corresponds to the classical separatrix. 
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FIG. 4. Time dependence of the energy dispersion for 
H — 1.25 ■ 10~ 4 and initial states: (1) near the center of cou- 
pling resonance, (2) and (3) - above and on the separatrix. 



corresponds to the Shuryak border |9(], the latter that 
only one quantum state is inside the classical stochastic 
layer. In this case, the classical chaos is completely sup- 
pressed by quantum effects. We have also found that for 
large times, t > 1000T, the diffusion-like growth of the 
energy terminates, thus indicating the localization of the 
classical diffusion along stochastic layers. 

In conclusion, we have shown that when the number 
M s of stationary states inside classical chaotic separatrix 
layers is relatively large, one can observe weak quantum 
diffusion that is similar to the classical Arnol'd diffusion. 
In our model this diffusion occurs along the nonlinear 
coupling resonance, due to the influence of guiding reso- 
nances originated from the external time-dependent per- 
turbation. Our results may find a confirmation in experi- 
ments on one-electron dynamics in 2D quantum billiards 
with time-periodic electric fields. 

This work was supported by grants RFBR No. 01- 
02-17102 and Ministry of Education of Russian Feder- 
ation No. E00-3. 1-413. FMI acknowledges the support 
by CONACyT (Mexico) Grant No. 34668-E. 



We calculated quantum and classical diffusion coeffi- 
cients and found that the diffusion in the quantum model 
is systematically weaker than the classical Arnol'd diffu- 
sion, see Fig.||. However, the global dependence of the 
quantum diffusion coefficient on the control parameter 
1/ yfjl is similar to the classical one. The data in this 
figure are obtained for initial states corresponding to the 
centre of the stochastic layer, where the diffusion coeffi- 
cient is maximal. For initial states corresponding to the 
border of the stochastic layer, the diffusion coefficient 
strongly fluctuates and can be several times smaller. 
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FIG. 5. Quantum (squares) and classical (solid line) diffu- 
sion coefficient vs the parameter 

To count the number M s of eigenstates inside the clas- 
sical separatrix layer, we calculated the width of the 
classical layer numerically, together with the number of 
energy levels in the corresponding energy interval. For 
fi > 1.25 • 10~ 4 the number M s was found to be about 
10. The right square in Fig.|5| for the lowest value of fi 
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